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DISCUSSION OF THE NON-ELASTIC BEHAVIOR 
OF BRIDGES UNDER IMPULSIVE LOADS 
PROCEEDINGS-SEPARATE NO. 185 


Rosert J. Hansen,'® A. M. ASCE, and Joun S. Arcuer,'* J. M. ASCE.— 


'S Associate Prof, of Structural Eng., Dept. of Civ. and San. Eng., Massachusetts Inst. of Technology, 
Cambridge, Mass. 

6 Asst. Prof. of Structural Eng., Dept. of Civ. and San. Eng., Massachusetts Inst. of Technology, 
Cambridge, Mass. 


The authors’ paper is of considerable value to engineers concerned with the 
vulnerability of bridges to attack from atomic weapons. The analysis as 
presented is limited to the case of attack from directly above the bridge and 
only mentions the possibility of attack from weapons that are detonated to 
either side or end of the bridge. For a thorough assessment of vulnerability, 
a complete consideration should be made of attack situations in which the blast 
wave approaches the bridge from the side. In many cases this situation would 
be most likely to cause failure of the bridge. The procedures of analysis 
presented by the authors, however, can be utilized for the latter conditions. 

Several factors not mentioned by Messrs. Fraenkel and Grinter deserve 
consideration. First—as regards the side loading condition—the possibility 
of sliding a bridge off its piers through failure of the anchor bolts is a mechanism 
of failure which deserves investigation. Second, truss members in compression 
that are subject to lateral loading from the air blast are particularly vulnerable 
if the axis of the member is normal to the direction of the air blast. The lateral- 
blast loading enhances the possibility of buckling of the member under its 
axial load. 

It should be apparent that the blast-load impulse is dependent on the 
percentage of openings in the bridge deck as well as on the width of the deck. 
The load impulse transmitted to an open-deck (steel-grid) highway bridge or 
an open-deck railroad bridge may be so small as to leave the bridge undamaged 
from a nominal atomic bomb burst 2,000 ft in the air over the bridge, without 
regard to the height of the bridge above the ground. 

The engineer should have a thorough background"’:'*."*..2t on the subject 


Shock Loading of Rectangular Structures,” by Walter Bleakney, Technical Report II-11.  Prince- 
ton University Press, Princeton, N. J., January, 1952. 

"4 Simple Method for Evaluating Blast Effects on Buildings," Armour Research Foundation of 
Illinois Inst. of Technology, Technology Center, Chicago, Hl., July, 1952. 

“Wind Loads on Truss Bridges,” by J. M. Biggs, Proceedings-Separate No. 201, ASCE, July, 1953. 

*~ ‘Proceedings of the Conference on Building in the Atomic Age,"’ Dept. of Civ. and San. Eng., 
Massachusetts Inst. of Technology, Cambridge, Mass., June, 1952. 

“Proceedings of the Symposium of Earthquake and Blast Effects on Structures,” Earthquake Eng. 
Research Inst. and Univ. of California, Los Angeles, Calif., June, 1952. 


of air-blast loading and dynamic analysis before attempting to analyze bridges 
for the type of behavior mentioned by the authors. 


8. J. FRAENKEL” anp L. E. Grinter,* M. ASCE.—The interest of Messrs. 


* Mer. Propulsion and Structural Research Dept., Armour Research Foundation of Illinois {nst. of 
Tech., Chicago, Ill. 
23 Dean of Graduate Scheol, and Director of Research, Univ. of Florida, Gainesville, Fla. 


Hansen and Archer in bridges subjected to impulsive loads is greatly appreci- 
ated. The writers are in agreement with the points cited by the discussers. 

It is entirely possible that atomic explosions resulting primarily in the lateral 
loading of bridges may have greater military effect than explosions which 
occur directly over a bridge. Lateral loading would be much more likely, for 
example, to produce alternate types of damage to the bridge—such as rigid- 
body tumbling or the sliding of the bridge off the supporting pier. As it was 
not the writers’ intent to establish a military doctrine for attack on bridges, 
consideration of these other possible causes of damage to bridge structures 
was purposely ignored. 

It is evident that openings in the bridge deck will have an effect on the 
damage sustained. Quantitative extrapolation to open-deck bridges cannot 
be made from the material presented. However, there is reason to believe 
that the general relations derived (in terms of bridge length, width, and height 
above ground) will also apply to open-deck bridges. 

Two items are worthy of additional emphasis. The first is the effect of 
longitudinal inertia—which is included in the equations characterizing bridge 
behavior, and which is the principal feature distinguishing this paper from the 
usual single-degree-of-freedom treatment of structural vulnerability. The 
importance of the longitudinal inertia is approximately the same as that of the 
rotary inertia, and it does not seem generally permissible to omit the longi- 
tudinal-inertia term from consideration. The second item relates to the failure 
criterion (the ratio of deflection to span length). The criterion is admittedly 
simple, but it does have the great advantage of lending itself to computation 
and of having significance (to a greater or lesser degree) for all types of bridges. 
It is possible that for given bridge types, and for given functional purposes, 
more pertinent and more specific criteria could be given. If this were done, 
the establishment of general rules of behavior would be rendered most difficult 
(if not impossible), and not even an approximate generalized treatment could 
be given to this subject. 

It should be realized that it was the writers’ purpose to develop and illus- 
trate a convenient method of handling the bridge-vulnerability problem and 
to derive certain general relations applicable to all bridges. It was not their 
intent to provide quantitative information on the vulnerability of any specific 
bridge. As a result of the continually changing state of knowledge con- 
cerning blast loading, quantitative answers for specific bridges should always 
be obtained with reference to the most recent information. 
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DISCUSSION OF APPLICATIONS OF THE RELAXATION TECHNIQUE 
IN FLUID MECHANICS 
PROCEEDINGS-SEPARATE NO. 223 


MLADEN BORELE®*—As indicated by the variety of results contained in the 
author’s paper, the relaxation method is a powerful means for solving engin- 
eering problems. Since the method is essentially one of successive approxi- 
mation, so that the process involved is of necessity rather laborious, it is of 
great importance to explore all possibilities for reducing the necessary time 
to attain a certain specified accuracy. From the experience gained in the cal- 
culation for free-surface flow toward partially penetrating wells by the relax- 
ation method, which the writer undertook at the Polytechnic Institute of the 
University of Grenoble under the direction of Professor Hunter Rouse, he be- 
lieves that the fineness of the mesh and hence the labor involved depend to a 
great degree upon the treatment of the irregular stars and the singular points. 

There are three ways in which irregular stars can occur. They appear (A) 
in the region where a finer mesh matches with a coarser one or (B) at a fixed 
or free boundary, or else (C) they can simply be chosen under favorable con- 
ditions as a matter of preference. Wherever irregular stars are encountered, 
it is necessary to have a formula for the calculation of residues. For an ir- 
regular star having legs of different lengths as shown in Fig. 1, one can obtain 
from Eq. (7) of the authors’ paper the formula 


for the residue, whereas, using Taylor’s expansion, the writer obtains 


g)+ + (2) 


which, in spite of its apparently more complicated form, is no more difficult 
to apply than Eq. (1). A comparison of these two formulas in connection with 
each of the three ways in which the irregular stars can occur will show their 
relative merits. 
A. Transition from a coarser to a finer mesh 

For axisymmetric flow, the formula for the residue at the joint of a regular 
star is 


a. Assistant, Fluid Mechanics Laboratory, Polytechnic Institute, University of 
Grenoble, France. On Leave, Hydrotechnic Institute, Serbian Academy of 


Sciences, Beograd, Yougoslavia. 
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in which the authors’ notations are used. The quantities at the right of Eq. (3) 
involving the derivatives of y will be denoted by AE. It is immediately clear 
that if these derivatives are large and r is small, AE is not necessarily 
small. Thus the greater the variation of the function and the smaller the 
values of r, the finer the mesh must be. This seems to have been universally 
recognized, for previous works have shown that the fineness of the actual mesh 
used varies from region to region. Since the transition from a fine mesh to a 
coarser one is effected by irregular stars, it is evident that if the calculation 
for the irregular stars is not very precise it is necessary to use a finer mesh, 
thus increasing the number of joints. 

The case of axisymmetric flow into a confined well, for which the exact 
solution is known, can be used as a basis for comparison. Supposing that 
>,- 1 and 9ns1 in Fig. 2 are the correct values, one can calculate $, accord- 


ing to Eqs. 1 and 2. From %,_, and $,,1 with the value of ¢, obtained in 
each case, it is possible to wanes the discharge per unit length q, through 
a cylindrical — of radius (n - pa and q, through a cylindrical surface 
of radius (n + 5 1). These stein in general differ from the correct value 


q. The one, 2 values of the deviations of q, and q, from q as calculated 
from Eq. (1) are denoted by 4q, and 4q,, and their mean value is denoted by 
4q°* The corresponding quantity calculated from Eq. (2) is denoted by Aq’. 

The following table shows that for a specified accuracy the value of n corre- 
sponding to Eq. (1) is larger than that corresponding to Eq. (2), so that Eq. (2) 
is preferable to Eq. (1). 


B. Boundary conditions 
Irregular stars usually occur at the boundary, which may be fixed or free. 
For a free surface the boundary condition can only be checked after all the 
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residues have been annulled. If the expression for the residue is not correct, 
the free-surface position will not be correctly determined. 

To make a further comparison of Eqs. (1) and (2),one may consider the 
region near the free surface in seepage flow. Beside satisfying the Laplace 
equation, the potential function ¢ has to satisfy the following two conditions at 


the free surface: (a) ? = z and (b) ge = 0, in which z is the elevation of the 


free surface, and n is measured in a direction normal to the free surface. Be- 
cause of (b) it is possible to consider the point A’ instead of the point A (Fig. 
3), thus eliminating one irregular leg and leaving only one leg of irregular 
length. Equation (1) gives 


2 
which becomes 


a + (5) 


Similarly, Eq. (2) gives 


2 
B= 


2 2 
(6) 


To make a cursory comparison of Eqs. (5) and (6), one can let 


A,;—>O ) sin xX —+>O 


In the limit, Eqs. (5) and (6) should give directly the formula for the residue 
at a joint located at a horizontal impermeable bed. This is the case with 
Eq. (6), but not with Eq. (5). The coefficient 2 appearing in the left side of 
(6) shows that the direction normal to the free surface is a privileged one— 
a property which is not displayed by Eq. (5). Since Eq. (5) comes from Eq. (1) 
and Eq. (6) from Eq. (2), this means that, once again, Eq. (2) is preferable to 
Eq. (1). 

The method of reflection used in Fig. 3 is not always convenient. One may 
use the ratio 2:1 for the lengths of the legs. If the joint is located on the sur- 
face, Eq. (2) gives 
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whereas Eq. (1) gives 
7 
a V0 $+ P-2 $+ a astra + (8) 


Since the numerators of the fractions are approximately equal, the residue 
calculated from Eq. (7) is twice that calculated from Eq. (8). 


C. Rectangular meshes 

In case the gradient of the unknown function can be judged a priori to be 
considerably larger in one direction than in the other, rectangular meshes 
can be used advantageously, with the long legs in the direction along which the 
gradient is small. Sometimes the number of joints involved can thus be re- 
duced by 30 to 40 percent. Evidently a more exact formula for the residue 
will have here a considerable importance. 

If, under favorable conditions, the mesh is made rectangular, with the longer 
- twice as long as the shorter ones, Eq. (2) gives 


whereas the left side of the corresponding formula given by Eq. (1) is only half 
as great. The difference is, consequently, very large. 

After making the foregoing comments on the formulas for the residue at the 
joint of an irregular star, the writer wishes to discuss a method of treating the 
singular points in some detail. The method involved consists essentially in the 
proper weighting of the legs, with the purpose of avoiding an extremely fine 
mesh which would otherwise be necessary near the singular points. 

In Fig. 4, after the values of the potential are known at the points around 
A, B, and C, one can form the residues, which will be seen to differ consider- 
ably from zero. Since these differences arise largely from the legs containing 
the singularity, one may choose certain weights for these legs so that the resi- 
due vanishes. For the point B the difference ¢, - ¢p must be multiplied by 
0.77 instead of 1. The weights for the legs of the various joints are shown in 
Fig. 4. 

To give an idea how the method of weighting can be applied to the calculation 
for a well, one considers the potential function 


p= bogr+N 


One has, in this case, 
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in which 


k= ( (2)- m= Bate (/2) 
For the case shown in Fig. 2, one obtains the coefficients presented in 
Table 2. 
TABLE 2 


In conclusion, the writer wishes to express the opinion that in each of the 
three cases in which irregular stars can occur Eq. (2) is preferable to Eq. (1). 
Moreover, by the use of proper weights considerable economy in labor can 
also be achieved in the vicinity of a singular point. 


2 
Ay 1.7778 1.5788 1.5031 1.4628 

| a7 | 2205s | | 

396-7 


| 
| et 
f, 
| 
| b 
| | 
| ¢ 
R=na | 
Figure 1 | 
| | 
| a 
| a | 
R=na 
| 
| 
| 
Figure 2 
396-8 


Figure 3 


Figure 4 


\ \ 

| A Q Z= 

| 
7 
B 

(I) (I) 
Z (2) (077 
B’ A Sc 
(2) (2) (1) 
(O77) 
(1) (1) 
Cc 
(1) 
396-9 


TURGUT SARPKAYA? J.M. ASCE—The engineering profession is indebted 
to the authors for their informative presentation of the relaxation technique as 
it is applied to problems in fluid mechanics. Numerous examples of a practi- 
cal nature are a valuable feature of the paper. Having worked under the super- 
vision of the authors in the solution of several such problems, the writer wishes 
to emphasize that both the direct and reverse relaxation methods, although by 
nature approximate and relatively time-*:onsuming, have distinct advantages 
over such approximate methods as Von Karman’s source-sink method and the 
flow-net method. In addition, the writer would like to enlarge slightly on their 
presentation by commenting on (1) the applicability of the reverse method, 

(2) the application of the direct method to have problems, (3) the accuracy of 
results, and (4) the treatment of sharp corners and singularities. 

As was emphasized in the paper, the application of the direct relaxation 
method, using the x-y plane, to the problems in which some part of the boun- 
dary is exposed to constant pressure is not a rational approach to the answer, 
since it is very difficult to estimate to what extent the boundary should be 
changed between successive trials. However, the application of the reverse 
method, especially to free-jet problems, does provide a rational approach to 
the solution. The reverse method, as it is presented in the paper, has the ad- 
vantage that the assumption of the variation of the potential with distance, i.e., 
the assumption of an x(¢)-curve, is made only once. The relaxation procedure 
for y is then completed, the results are checked against one of the Cauchy- 
Riemann equations, and a new x(?%)-curve is obtained. This procedure elimi- 
nates the necessity of assuming another entirely arbitrary x(%) relationship. 
The writer’s experience leads him to recommend that the x-values for the 
final trial be computed with a relaxation similar to that for the y-values, 
rather than determined directly from Cauchy-Riemann conditions. 

The authors point out that the relaxation method is applicable to certain types 
vf wave motion. An example, the case of irrotational standing waves in a chan- 
nel of constant depth and variable width, will be discussed after the following 
explanation of the general procedure. Many types of wave motion are essen- 
tially irrotational, and the motion can therefore be defined in terms of a veloc- 
ity potential. 


p= cos cash k (h+Z) (1) 


in which ¢ satisfies the Laplacian equation, 
2 
(2) 
Also, if the amplitude of the motion is very small, the well-known condition for 


the free surface is 
o¢ a = 


in which z is taken as positive in the upward vertical direction. If the value of 
? from Eq. (1) is substituted into Eq. (2), / (x,y) is found to satisfy the cri- 


terion, 


b. Research Assistant, Iowa Institute of Hydraulic Research, Iowa City, Iowa 
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Also, if Eq. (3) is to be satisfied, 
2 
an 
(4 k 9g tanh kh =O (5) 


The two relationships, Eqs. (4) and (5), are fundamental to this study as well 
as to studies of other physical problems. The goal of the analysis is the at- 
tainment of a solution that satisfies Eq. (4) and the boundary condition, 
09/On = O. Equation (5) represents the relationship between the period of 
the wave, the depth of the water, and the parameter k. Equation (4) can be 
written in the usual finite-difference form, 


_— with k taken from Eq. (5) by trial, arbitrary yw values assigned, and 


values computed, the relaxation process can .be performed. After it 


i=1 
is completed, the elevation " of the free surface is obtained from 


The troughs and crests of succeeding waves are obtained by connecting, re- 
spectively, locally maximum and minimum 1 -values. 

An application of this result is presented in Fig. 1. A parabolic transition 
reduces the width of a rectangular, constant-depth channel to one fourth of its 
original width, as shown. The depth of the still water in the channel was as- 
sumed to be 0.6 ft, k = 7 ft™!, and the period 0.64 sec. Figure 1 shows the 
instantaneous pattern of the crests of standing waves in the channel. The am- 
plitude of the wave, also indicated on the figure, was found to be inversely 
proportional to the square root of the width, in accordance with Green’s law 
[10]. If constant depth can be assumed, the relaxation method can be applied 
in the analysis of seiche in irregularly shaped harbors, for which standard 
methods of analysis fail. 

The authors caution that the accuracy of a given result can only be judged 
by computing the velocity or pressure distribution along a boundary or other 
arbitrarily selected lines. However, the computation of velocities on a free 
boundary even from a fine mesh if the mesh size near the boundary is of the 
order of magnitude of the radius of curvature of the boundary will not lead to 
an absolute check because of the variation of A-values in the star. If a star 
with a very small A-value is to be considered, a formula for curvalinear ex- 
trapolation must be used. 

The solution of certain practical problems may present very difficult boun- 
dary conditions, such as sharp corners which give rise to singularities of 
various types. Special treatment is essential if the function approaches infinite 
values near the corner. This difficulty can be overcome by the application ofa 
somewhat more complicated metinod. 

The solution of V?%=0O can be expressed by the series 
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in which n = a (k = O, + 1, + 2, ...) and @ is the angle of the sharp corner 
as is shown in Fig. 2. For & = 271, @ takes the form 
32. 3 
p=A,+A,r +A, rcose+A,r O+ (9) 


Only the first four terms need be retained. The units of r can be so chosen 
that r=1 at A, B,C,and D. Intermsof ¢,, ?p, %c, and $p, which 


are known, one obtains the coefficients A,, A,, A,, and A,. Using these co- 
efficients, Pat B1° $c 1, and %p can be expressed in terms of ‘ 


and p, as follows: 
4572 g +.23519 +.20936 +.0992 
.235/ g + .5930 .0992 +.0728g 


=. /O 
2093 + .0992 $ +.4572$ + .23519 (/0) 


=.0992 ¥ + .0728 +.235/ g + - 5930 


The mathematical analysis of the above method can be found in references 
[11] and [12]. These are the equations used at singular points. The relaxa- 
tion procedure is carried out in the usual manner, observing that Eqs. (10) 
hold at singular points. 

The usefulness of the relaxation method in a wide variety of problems, as 
well as for a variety of specific boundary conditions for a specific problem, 
has been demonstrated. 
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DISCUSSION OF NUMERICAL ANALYSIS 
CF FRAMES WITH CURVED GIRDERS 
PROCEEDINGS--SEPARATE NO. 250 


WILLIAM A. CONWELL’, M. ASCE—A wealth of data which will be of great 
value to structural engineers is included in this paper. It is the purpose of 
this discussion, after some general comments, to indicate points which will 
demonstrate the wide application of the informatior given in the paper. 

In the first paragraph under the heading “Stiffness and Carry-Over Factors” 
the author does not state that, during the application of the moments, the ends 
of the curved girders are restrained against translation although he obviously 
intends this. It is believed that the paper should be explicit on this point par- 
ticularly in view of the current work? by L. E. Grinter, M. ASCE and C. H. Tsao 
which deals with moment distribution where translation of the joints is per- 
mitted. 

Though the author has devised excellent nomenclature for his work, the 
writer would prefer to see a slightly different expression for the carry-over 
factor. The author chooses to use the first sub-script inside the parentheses 
as indicating the axis to which the moment is carried; the second sub-script 
to indicate that about which the moment is applied. Interchanging these two 
sub-scripts would put the sub-scripts inside the parentheses in parallel with 
those outside. For instance, (Cxz)p_a would then be the carry-over factor 
from the X axis at.B to the Z axis at A. The author, no doubt, has very good 
reasons for using his nomenclature. It would be of value to know what lead to 
its adoption. 

The curves given in Figure 3 are most enlightening and are a fine repre- 
sentation of tabular data. The additional curves, to which the author refers 
under the heading “Stiffness and Carry-Over Factors,” may have sufficient 
value to be included in the closing discussion. 

One of the beauties of the paper is the fact that it permits the same proce- 
dures to be followed for circular girders that we have been in the habit of using 
for straight beams.° For instance, if one desires to obtain the fixed end mo- 
ment for a uniformly distributed load on a part of the span of a circular girder 
one may initially solve the continuous girder shown in Figure 9(a) with a false 
support at C starting with fixed-end moments of Table 4. After computing the 
reaction at C one can, from Table 3, obtain the fixed end moments for a load 
W at C equal to and opposite in direction from the reaction (see Figure 9(b)). 
The super-position of Figures $(a) and 9(b) yields Figure 9(c) which shows a 
circular girder partially loaded. It is recognized that Table 3 of the paper can 
be used to obtain this same information by taking a series of positions for the 


1. General Engineer, Structural Engineering and Design Department, Duquesne 
Light Company, Pittsburgh, Pennsylvania. 

2. Separate No. 298—Proceedings of American Society of Civil Engineers, 
October, 1953, entitled “Joint Translation by Cantilever Moment Distribu- 
tion” by L. E. Grinter, M. ASCE and C. H. Tsao. 

. “Continuous Frames of Reinforced Concrete” by Hardy Cross and Newlin 
Dolbey Morgan, John Wiley and Sons, Inc., New York, N. Y. 1932. 
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load W. On occasions when, as stated under “Determination of Stiffness, 
Carry-Over Factors and Fixed End Moments,” it is desired to consider shapes 
of girders other than circular for which data similar to the author’s is not 
available, the procedure just outlined may result in a saving of time and effort. 

When one thinks of the over-all problems of Numerical Analysis of Frames, 
there immediately comes to mind not only the fixed end moments, stiffness 
and carry-over factors mentioned by the author, but the effect of translations 
of the joints. In the treatment of circular girders the problem could enter in 
analyzing the effect of settlement at a column or perhaps, in extreme cases, 
for providing for the axial deformation which might take place in a compara- 
tively light column or hanger. Although the author does not specifically men- 
tion translations of joints he has provided the data for computing this effect. 

As with straight members it is apparent that the effect of translation of the 
joints can be obtained by introducing equal rotations at each end of a member. 
For instance, as shown in Figure 10, if end A of a circular girder were to set- 
tle an amount A the effect could be duplicated by adding the effects of rota- 
tions around Z axis introduced at ends A and B and equal to 4/r. Figure 10 (a) 
shows the effect of rotating end A through an angle of 4/r and Figure 10 (b) 
shows the same for an equal rotation of end B. The superposition of Fig- 
ures 10(a) and 10(b) gives Figure 10(c) which shows the fixed end moments 
due to a settlement or translation of A at end A, while A is held against rota- 
tion and B against both translation and rotation. These fixed end moments 
may be distributed as any other fixed end moments to determine the over-all 
effect on a structure of which member AB is a part. 

The writer is pleased to have had the opportunity to discuss this paper, 
which he considers to be of high calibre and of great value to the engineering 
profession. 


Corrections to Proceedings-Separate No. 250. 


In Table 2, the value of the carry-over factor for EI/GJ=1, ¢ = 60°, 
and 7 =O should be, inthe column headed Kz, (Cxz) _ 4 (t/EI) 


= Kya (Cza)a-a(r/ED, equal to -0.482. In Table 5, the value of 
(Cy x)3 Should be -1.000. In Table 6, the value of (Czz)3_9 should 
be +0.487. 
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DISCUSSION OF MECHANICS OF MANIFOLD FLOW 
PROCEEDINGS-SEPARATE NO. 258 


JAMES R. BARTON.—An excellent job has been done by Dr. McNown in 
summarizing the basic problem of manifolds in a clear and logical manner, 
and engineers working on the design of projects involving manifold flow 
should benefit from the information presented. In a paper of such wide scope, 
one finds it difficult to present all of the details connected with the research 
even though some of them may be pertinent to an evaluation of the data. 
Therefore, it may be profitable at this time to briefly describe some of the 
experimental techniques and to evaluate the accuracy of the measurements 
for +: one work on dividing flow as performed by the writer in 
1946 [1]. 

All discharge measurements were made in weighing tanks and a minimum 
of two measurements were made for each discharge. A careful study of these 
data indicated a maximum probable error of about 1 percent. However, the 
measurements of piezometric head were not quite as accurate. Although cap- 
illary tubes were inserted in the lines to damp out fluctuations in both the 
mercury manometer and the air-water differential gage, some fluctuation 
still existed. A study of the piezometric head data resulted in an estimated 
maximum probable error of about 2 per cent with a possibility of 3 per cent 
error in the vicinity of the lateral where excessive disturbance occurred. 

Determination of the proper hydraulic gradient for each pipe was at first 
a problem, since the individual piezometric head measurements showed some 
scatter. In order to bypass the repeating decision of “what should the exact 
slope be on this run?”, each lateral and the main pipe were carefuily cali- 
brated by making a log plot of Reynolds number vs. hydraulic gradient. Each 
of these plots was based on a minimum of nine different slopes which in turn 
were based on a minimum of 6 point measurements of piezometric head. 
With these calibration plots the correct slope could always be fitted to the 
piezometric head data for individual runs even though the data exhibited some 
scatter. Of course, in the case of excessive scatter the run was discarded 
and repeated. A study of Fig. 3 will indicate that a change of slope of any of 
the hydraulic grade lines could affect the terms Ah and hz, especially when 
these terms were small. For this reason the precaution of using the correct 
slope for each flow measurement is believed to have avoided some errors 
which might have resulted in a larger scatter of the data plotted in Figs. 4 
and 5. 

In order to illustrate the usefulness of the author’s data as a guide in re- 
search or practice, the writer would like to briefly present some findings of 
a Master’s research project completed at the University of Utah in 1949 by 
Mr. Grant Borg.! With the object of establishing a workable design for the 
underdrainage system of small rapid sand filters, Mr. Borg ran extensive 
tests on the discharge characteristics of a perforated pipe underdrain. The 
research resulted in the establishment of certain relationships which must 


1. “Design Data on Perforated-Pipe Underdrains for Small Rapid Sand Fil- 
ters” by Grant Kenneth Borg, M.S. Thesis, University of Utah, 1949. 
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occur in order to give satisfactory distribution of wash water through the 
perforated pipe system. 

The equipment consisted of one-4-inch manifold 6 feet long with 12 one- 
inch laterals attached to each side of the pipe on 6-inch centers so that all 
laterals were in a single plane. Orifices varying from 1/8 inch diameter to 
7/16 inch diameter were drilled in the bottom of the laterals on 4.18 inch 
centers. In this way, the 4-inch pipe acted as a manifold for the lateral pipes 
which extended about 21 inches from the. center of the 4-inch pipe, and each 
lateral in turn was a manifold system for the orifices which discharged into 
the:atmosphere. Since the orifice discharge was not always at an angle of 
90° to the center line of the lateral, Mr. Borg’s data cannot be compared in 
the strictest sense with the author’s data, but in some details a qualitative 
comparison exists which is interesting to note. 

Although other comparisons might be made, only one will be included here 
to illustrate how the author’s data would have been useful in the research at 


Utah. For small a ratios, and varying @ ratios the flow distribution 


through the orifices was very uniform as Fig. 5 indicates it should be. For 


” = 1/4, the ratio of a remains essentially constant over a wide 
V; /2g 
1 


range of discharge ratios Q . However, as 5. increases, the ratio of 


Q D 
hi 
: Q 

7s, varies considerably over a wide range of—+. The results of Utah 
Q 
clearly illustrated that as the orifice was enlarged, the flow distribution be- 

hf 
Vj /28 


and Q,/Q naturally varied greatly from the first orifice on the lateral to the 
last one. The laboratory research at Utah was for a specific purpose, and 
the work was very well done, but the author’s results could have been effec- 
tively employed as a guide, thereby reducing the total amount of work re- 
quired for a significant answer. 


came poorer because the ratio of varied considerably with Q,/Q 


| | 
| 
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M. R. CARSTENS?! and JOHN ROBERSON?— Professor McNown has clearly 
and concisely described the mechanics of manifold flow. The results present- 
ed in this paper are a demonstration of the values to the engineering profes- 
sion of wisely directed graduate-student research programs. 

With the principal characteristics of the flow defined for a single branch 
point in a conduit, it is possible to incorporate the results into the analysis of 
a multiple branching conduit system. The only additional restriction ‘to such 
an analysis is that the branch points be separated by a sufficient length of 
straight conduit. The reported experiments are indicative that a spacing of 
20 diameters between branch points is sufficient. As mentioned by Professor 
McNown, the multiple branch analysis must be started with the farthest down- 
stream branch and must be continued by successively accounting for the ef- 
fects of each section of intervening main line conduit and of each branch. The 


Figure 1. Definition Sketch for Multiple Branch Flow. 


general expression relating the piezometric head in the upstream or n th 
branch to the downstream or 1 st branch (figure 1 and figures 1a and 3a, 
McNown) is as follows: 


ah 
Wey}, 29 


h-h, Va 
J, 
VW 
J, *9 


(1) 


1. Associate Professor, School of Civil Engineering Georgia Institute of Tech- 


nology, Atlanta, Georgia. 
2. Assistant Hydraulic Engineer, Division of Industrial Research and Instruc- 
tor, Civil Engineering, State College of Washington, Pullman, Washington. 
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The parameter fL/D is the loss coefficient for the straight pipe section 


between the branch points. The values of 4 h/V?/2g can be determined direct- 


ly from figure 4, McNown. The value of (h -h,)/V?/2g can be determined 
from figure 5, McNown, as follows: 


Vg * B) (2) 


The computation procedure indicated by equation (1) is satisfactory for an- 
alysis of an installed system but laborious trial computations are involved in 
the calculation for design. 

In many situations the multiple branch variables are limited and the solu- 
tion of equation (1) can be presented on a nomograph. Work was done at the 
State College of Washington, Division of Industrial Research, by the writers 


to develop the nomograph of equation (1). The multiple branch simplifications, 


facilitating a graphical solution, are chosen as follows: (a) all branch points 


are identical; (b) all sections of main line conduit connecting the branch points 


are identical; and (c) all branch lines are identical both as to diameter D, and 
as to coefficient of discharge Cg in the expression 


With these simplifications equation (1) can be reduced as follows: 


+ {Mey * (4) 


Mj is the ratio of the discharge of the ith branch to that of the downstream 
branch, i.e., Mj; = Q;/Q). Thus My is the ratio of the discharge in the up- 


stream branch to the discharge of the downstream branch. 
Equation (4) has been solved for various values of the dimensionless 
parameters Cy, {L/D, D/ D), and n. Figures 2 and 3 are graphical presen- 


The value £ M/n is the ratio of the average discharge of all branches to the 
discharge of the downstream branch, that is, Qaye / Q;. The nomograph is 
used by following horizontally or vertically from curve to curve of known 


(3) 


and 
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Figure 3. Graphical Solution for Multiple Branch Flow. 


value until the scale or curve of unknown value is reached. This procedure is 
illustrated by the light dashed curve on figure 2 from which the values M= 1.58 
and 2 M/n = 1.16 are determined for n = 25, D/D, = 3.0, {L/D =6 and 

Cq = 0.06. 

The graphical solution (figure 3) can be used for the design of portable 
sprinkler irrigation lines. The portable sprinkler irrigation line consists of 
a series of light metal tubes which are inserted into quick-coupling connectors. 
A riser tube (generally 1-inch) protrudes from the top of the connector. A 
rotating sprinkler is mounted on top of the connector. Figure 4 is a photo- 
graph of a typical line. In order to use figure 3 for sprinkler irrigation line 
design, there are certain additional approximations and restrictions that will 
be discussed in the following paragraphs. 

The first approximation concerns the proper choice of the constant {L/D 
curve on figure 3. The tubing sections of a portable irrigation line are gen- 
erally of extruded aluminum. As the roughness height is very small, the new 
tubing sections would be classified as “hydraulically smooth.” Hence the 
value of the Darcy-Weisbach friction factor f will be a function of the Rey- 
nolds number of each tubing section. The Reynolds number will be the least 
value in the downstream section and will increase successively to the 
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Figure 4. Portable Sprinkler Irrigation Line. 


greatest value in the upstream section. The best value of a constant f to re- 
place the variable f can be determined by a summation of the (fL/D) (V2/2g) 
terms. From this summation, the average velocity VR to be used in calculat- 
ing the Reynolds number is found to be 


yy (5) 


Thus, one can determine the best value of the Reynolds number by using VR 
By referring to standard curve of f versus Reynold’s number, the best vale 
of f can be determined. The parameter {L/D is then computed using this 
best value of f along with the properties of the repeating tubing sections L 
and D. 

In sprinkler irrigation line design the use of the sprinkler coefficient of 
discharge is much simpler than the use of the branch line coefficient dis- 
charge. In equation (1) the difference hjp - hl, was expressed in terms of 
flow and geometric properties. In equation (4) the values of hj, and hj, 
were eliminated, as previously explained, by the use of the branch discharge 


equation, that is, 
2 2 
= (t ) A 
Cy/ 29 


Thus the value 


was used in the derivation of equation (4). However, in the case of sprinklers, 
the discharge is a function of the sprinkler geometry and the pressure at the 
base of the sprinkler 
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2 


in which the subscript s denotes sprinkler. It follows that 


+ + #2, 


(8) 


in which Zg is the elevation of the sprinkler base and the last term is the 
friction energy loss in the riser. The use of the sprinkler coefficient of dis- 
charge Cg, in place of the branch coefficient of discharge Cy is permissible 


if hj, - hj, = Psn/ ¥ - Ps,/¥. This condition is satisfied in a level irriga- 


tion line with equal-length risers and if the friction energy loss difference in 
the risers is negligible in magnitude. This difference is negligible for a num- 
ber of reasons but primarily because this energy loss is small in relation to 
the value of Pg/¥ used in irrigation lines. 

In summary, the results as reported by Professor McNown have been in- 
corporated into a graphical method for the solution of multiple branching con- 
duit systems. Nomographs (figure 2 and 3) have been prepared for systems 
in which, (a) all branch points are identical, (b) all sections of main line con- 
duit, connecting the branch points are identical, and (c) all branch linés are 
identical as to diameter D and coefficient of discharge Cy. Methods of 
utilizing the nomographs were discussed for the design of a level portable 
sprinkler irrigation line with equal-length risers. This discussion included 
the method of determining the best value of f to be used in the parameter 
{L/D and the feasibility of utilizing the sprinkler coefficient of discharge Cds 
in lieu of the branch line coefficient of discharge Cg. 

The results, as reported by Professor McNown, are deficient in respect 
to the minimum length of branch line. The analysis and experimental studies 
were predicated on the concept that the length of branch line was sufficient so 
as to have no influence upon the expansion of the jet (figure 2, McNown). It 
would seem possible to determine this minimum length from a study of the 
piezometric head line in the branch conduit utilizing the data obtained by 
Barton [1]. Professor McNown also fails to indicate why this data was not 
utilized in determining the minimum spacing between branch points in the 
main conduit. 


(7) 
2g 
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JULIO ESCOBAR,’ J. M. ASCE— The paper presented by Mr. McNown is 
representative of the kind of work conducted at the Iowa Institute of Hydraulic 
Research. Rather than consider particular problems from the standpoint of 
their immediate or practical applications, the Institute staff has been engaged 
in fundamental research into the mechanics of fluid motion. In evolving the 
different phases of the studies, new basic principles are giver. coherent ex- 
pression and their experimental foundaiions are thoroughly weighed and valued. 

Mr. McNown undertook such a basic investigation of the hydraulics of man- 
ifolds after the publication of the work of Soucek and Zelnick (11), together 
with his discussion thereof which merited the J. C. Stevens Award in 1946. 
After five years of research, Mr. McNown has published his findings. Faith- 
ful to the orientation of his Institute, he has presented not a series of graphs, 
tables, or charts applicable to the design of a particular type of manifold, lock, 
or water-supply system, but rather an interpretation of the behavior of flow 
in the vicinity of laterals, together with a definition of the general principles 
governing the hydraulics of manifolds. 

Among the many aspects that can be studied in the problem, the one per- 
taining to the explanation of the energy changes and the head losses that occur 
in the region of the lateral for dividing flow is very important. The writer 
was concerned with this phase of the investigation while at the University of 
Iowa, and wishes to contribute additional findings from his studies (3) which 
are believed to be helpful in the understanding of certain phases of the phe- 
nomenon. These subjects had been discussed briefly by other investigators 
(1, 5, 11), particularly in regard to the presence of negative losses in the re- 
gion of the lateral. 

The apparent increase in the energy of flow past the lateral for small dis- 
charge ratios (lateral to total) is explained by advancing the hypothesis that 
the lateral discharge comes from a region of low-velocity flow, containing 
below-average kinetic energy relative to the remaining flow. Therefore, the 
equation of energy, from which the equation 


was established by Mr. McNown in his paper, should be modified to include 
velocity correction coefficients? 
2 2 
oh + hy 


with reference to Fig. 1. 


lo, 


Figure 1 


. Assistant Engineer, R. J. Tipton y Asociados de Colombia Ltda., Bogota, 


Colombia. 
. Rouse, Hunter, Elementary Mechanics of Fluids. New York, John Wiley and 


Sons, Inc. 1946 -p. 114. 
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If the velocity correction factors are redefined as: 
Ky = I+€, and It 
and if Eq. (a) is rewritten in the form 


the pressure rise will be modified by two terms depending on the velocity 
heads and the velocity distribution upstream and downstream from the lateral 
if Ah, is the corrected pressure rise, 


ah, = 4h +6 (c) 


and the final equation takes the form, 


_ »_(%)?_ d 
(d) 


The second member of equation (d) is the same as that of McNown’s Eq. 
1 plotted in Fig. 4 of his paper. Therefore, the pressure rises may be 
computed and compared not only with the theoretical curves but with Barton’s 
experimental curves also plotted in Fig. 4- 

The remaining problem is the evaluation of the velocity correction factors, 
which are defined by the equation 


Ke= i (£)"dA (e) 


Their evaluation depends on the manner in which the main flow is divided in 
the region of the lateral. Once the section of the lateral discharge is known, 
the energy factors can be determined and applied to the general energy 
equation. 

It was necessary to observe in the laboratory what part of the cross section 
of the main pipe contributed flow to the lateral. A lucite “T” section, 2 inches 
by 1 inch, was inserted in the experimental 2-inch pipe to make possible ob- 
servation of the path followed by dyed droplets or by air bubbles injected at 
various points in a section well upstream. As the use of dye was somewhat 
troublesome because prolonged operation colored the recirculating supply and 
interfered with other experiments in the laboratory, most of the exploratory 
measurements were made with air. However, later analysis showed that the 
dye solution was more accurate for the intended purpose. 


Experiments showed that for values of Qu < 0.5 the area for Q was 
Q 


bounded by slightly curved lines convex to the side of the lateral as shown in 


Fig. 2. For Qi > 0.5 the experiments did not show exactly the form of the 
Q 


section but it seemed that the bounding curves increased their curvature to 
form almost closed boundaries as shown in Fig. 2. 
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Dividing lines 


Fig2-—LINES OF DEMARCATION BETWEEN LATERAL 
AND CONTINUING FLOW FOR THE FLOW 
APPROACHING THE LATERAL. 


In Fig. 3 are shown the results obtained from the tests with air and dye 
for a horizontal traverse of the hypodermic needle used for dye injection; the 
line shown was computed from the assumption, made for simplicity, that the 
area required for a given lateral discharge was limited by a vertical line as 
shown in the figure. This assumption was not exact but proved very helpful 
for the interpretation of results. Other traverses were made along lines 
inclined 30 and 60 to that for Fig. 3 as indicated in Fig. 2. 

It was mentioned before that air bubbles did not define precisely the boun- 
dary of the diverging flow as is borne out by the excessively high position of 
the points from the observation with air bubbles in the plot. This may be ex- 
plained from a consideration of the mass accelerations involved. In the re- 
gion of the lateral, the pressure gradient that produces the lateral discharge 
causes a larger acceleration of the injected air bubbles than of the surround- 
ing water because of the different densities. Although this relative accelera- 
tion is diminished by the resistance and virtual mass of the surrounding me- 
dium (the principle of the virtual mass is presented in most treatises on Hy- 
drodynamics**:5 the bubbles are carried toward the lateral pipe more rap- 
idly than the water or than the dyed droplets, the density of the latter being 
practically equal to that of the water. Consequently, the air bubbles do not 
follow the stream lines of the liquid, and, being displaced more than the fluid 
particles, the section of the lateral discharge determined from their paths 
appears to be larger than it should be. 

With the assumed form of the area for the portion of Q which is diverted 
to the lateral, and with the assumption of a velocity distribution described by 
the von Karman equation® 


3. Lamb, Horace, Hydrodynamics. New York: Dover Publications, 1945, p. 124. 
4. Green, S. L. Hydro and Aero-dynamics. London: Sir Isaac Pitman and Sons, 
Ltd. 1937, p. 114. 
5. Streeter, Victor L., Fluid Dynamics. New York: McGraw-Hill, 1948, p. 64. 
. Rouse, Hunter, Elementary Mechanics of Fluids. New York, John Wiley and 


Sons, Inc., 1948, p. 64. 
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Fig.3.— EXPERIMENTAL DETERMINATION OF THE 
LATERAL FLOW REGION-—For traverse 
angle of O°. 


= VE (2.15 log %, + 1.43) 4/ 


the energy factors were evaluated. From Eq. (d) it was then possible to com- 
pute corrected values of the pressure rise across the lateral, and compare 
computed values with those of Barton’s experiments. 


For > < 0.5 all recomputed points showed that the energy gain previous- 
ly found had been accounted for. For ae > 0.5 all recomputed points were 


abnormally high, indicating a failure in the assumption of the form of the area 
of Q, as has been explained before. However, some rise of the experimental 
petals is to be expected, because the water going to the lateral has above aver- 
age kinetic energy. Although it would have been possible to compute the cor- 
rect pressure rise following the exact flow pattern, it was thought that the 
refinements obtained would not have justified the work involved in the com- 


putations since high values of s are not common in practice. 


The apparent anomaly of the energy increase across the lateral for low 
values of Q) /Q was thus explained in a rational manner, supported by exper- 
imentation. The pressure rise is affected by the velocity distribution of the 
flow and either the theoretical curves of Mr. McNown’s Fig. 4 could be raised 
or the experimental curves lowered. 
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Fig.4.- COMPARISON OF OBSERVED AND COMPUTED 
EFFECTIVE CROSS-SECTION OF THE 
LATERAL . 


One second phase of the problem of lateral flow to be studied from knowl- 
edge of the flow pattern is that of the form losses occurring at the entrance 
to the lateral. To detect the eddy formed in the region of separation of the 
lateral flow, another technique was developed whereby dye was injected through 
a hole drilled into the lateral and the extent of the colored separation zone re- 
corded for various Q,/Q ratios. 

The reduction of the effective area of the lateral pipe due to separation of 
the flow as the water enters the lateral, is very similar to that occurring at 
an orifice plate or a flow nozzle. It seems logical therefore to compare the 
head losses in the lateral with those of the mentioned hydraulic transitions. 
Mr. McNown has shown the remarkable correspondence between the measured 
head losses and those based on theoretical considerations for plane irrotation- 
al flow. The writer also found good correspondence between his laboratory 
determination of the effective cross-section of the lateral and those computed 
by use of equations for the head loss for an abrupt enlargement (McNown’s 
Eq. (5)), and for a nozzle. For the latter, the empirical equation recommend- 
ed by Bean and Beitler’ 


hy = — (1-¢) 4) 


was used. This equation takes the dimensionless form: 


7. Bean, H. S. and Beitler, S. R., “Research on Flow Nozzles” Transactions 
ASME, Vol. 60, April 1938, p. 244. 
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Figure 4 shows the computed and observed values of the effective cross- 
section of the lateral for different Q,/ Q ratios. The contracted area decreas- 
es with Q,/Q for small values of the latter, thus increasing the head losses 
because of the greater deceleration experienced by the flow. For Q,/ Q = 0.35 
the area-ratio becomes sensibly constant and approximately equal to 0.5. 
These empirical findings are comparable both to the calculated and the ex- 
perimental values presented in the paper. 


JAMES J. HARTIGAN!—The theoretical relationship indicated by Eq. (4) 
was used in analyzing test data obtained by the Corps of Engineers in a model 
of a lock culvert with a single lateral port.? The coefficient of contraction as 
determined by the theoretical analysis should prove to be useful in the design 
of ports of lock manifolds. 

Head losses between the conduit and port were determined from the piezom- 
eter head plus velocity head in the conduit on the upstream side of the port and 
the chamber water surface plus the velocity head, Vp /2g, of the flow through 


the port based on the outlet area exclusive of the rounded edges. References 
to port entrance and outlet are defined by direction, in divided flow. Test data 
for combining flow are not included in this discussion. In a definition sketch 
such as Fig. 3a, the chamber water surface would replace the upstream end of 
the piezometric pressure gradient in the lateral. 

In Fig. A head losses in terms of V2, /2g are shown for four ports, three of 
which had rounded entrances and exits and the other square edges. Dimensions 
of the ports are shown in Fig. B. Although the port in each case was located 
at the bottom of the conduit, thus limiting the contraction in the vertical to the 
top side, there seemed to be no reason for changing the areas to be used in 
determining the ratio of a/b in Eq. (4). 

The ratio of port length to conduit width was evidently large enough so that 
values of C, computed for a barrier of infinite length were applicable. How- 
ever, it was observed in the model that some of the ports were too short to 
allow complete expansion of the jet before reaching the outlet. Test data are 
therefore compared with the theoretical curve for head loss in a sudden ex- 
pansion and also for the condition of no expansion of the jet, 


in which V} is the velocity of the jet. A curve midway between the ordinates 


of the curves for sudden expansion and no expansion is shown to help refer- 
ence the test points. 

The rounding of the entrance to the ports made indefinite the effective en- 
trance area for determining the ration a/b. Dimensions for port 2 in Fig. B 
show the width and height used in determining the maximum effective entrance 
area. In addition to the upstream and top curvature, the downstream curvature 
to the 45-degree radius is included. 

Since in three dimensional flow the contraction coefficient C, is applicable 


to an area, values of Vv C, will apply to the width and height of the port. For 
port 2 the jet would not be suppressed in the vertical until the value of C, 


3.24 
4.24 


1. Hydraulic Engineer, Corps of Engineers, St. Paul District. 
2. Hydraulic Laboratory Report No. 52, St. Paul District, Corps of Engineers, 
June 1946. 
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| . After determining the effective entrance area 
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for the conditions noted above, the value of Q,/Q = 0.7 was obtained from a 
graph of Eq. (4). For greater values of Q,,/ qd the port should flow full and 
head loss based on the contraction coefficient will be zero. For lower values 
of Q,,/Q the jet would be contracted in the vertical but not in plan until the 
maximum effective width can be used. 

According to the theoretical analysis the Wheeler port will have a contract- 
ec jet smaller than the minimum port area throughout the range of values of 
Q,/Q, thus the theoretical curves remain above the value for h'; = zero. 
Similar conditions exist for square-edged port No. 1. 

Data for port 5 are presented because the greatest range of the head loss 
ratio was measured in tests on this port. However, port 5 is probably the 
least practical of the four. It was assumed that the barrier in the plan view 
was of sufficient length to deflect the jet at an angle of 83 degrees. Accord- 
ingly a plot was made of Eq. (4) for the special case of the jet deflected 83 
degrees, However, for the range of available test data i.e., up to Q,/Q = 0.7 
it was found that values of C, were increased by less than 1 percent over cor- 
responding values for a 90-degree jet. Values of C, applicable to port 5 would 
probably vary from those for a 90-degree jet by less than 0.5 percent. 

Fig. A shows that only for port 5 do the head losses fall substantially above 
the midordinate line. For ports 1, 5, and Wheeler port there was some con- 
sistency in the relationship between the size of the jet and the outlet area of 
the port at the point where the head loss dropped below the midordinate curve. 

Flared ports 6 and 7 were more efficient than port 5. Flared 9 and 10.5 
degrees, respectively, these ports indicated head losses for values of Qp /Q 
greater than 0.15 which were not in excess of the loss that might be expected 
in gradual expansions twice as great. 


SADIQ M. NIAZ,! J. M. ASCE— As a part of the study outlined by the author, 
the writer made detailed observations on certain phases of converging mani- 
fold flow. In particular, information was obtained concerning the terms for 
the head loss and the unbalanced force for this type of flow. These observa- 
tions together with alternative forms of plotting are presented herein as an 

aid in the interpretation of the data. 

From application of a simplified energy relationship the author obtained 


Eq. (6): 


| 


Because the conduit velocity increases at the junction, the loss experienced 
by the conduit flow is somewhat similar to that for flow through a boundary 
contraction. It can, logically, be expressed in terms of loss coefficients to 
be applied either to the velocity head of approach or to the difference of the 

| upstream and downstream velocity heads; the two are simply related through 
| 


the discharge ratio and the diameter ratio. Thus, from the latter alternative, 
the head loss can be expressed in the form: 


2 2 
29 / yz (a) 


1. Director Central Designs, Punjab Irrigation Dept., Lahore, Pakistan. 
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in which m is a loss coefficient which varies with the aforementioned di- 
ameter and discharge ratios. 
Equation (6) can now be rewritten: 


A family of curves for various constant values of m relates the dimensionless 
pressure change and the discharge ratio. In Fig. A, representative curves for 
various values of m are compared with the data obtained for convergent flow, 
this plot being an alternative to the author’s Fig. 7. The curve for m = O cor- 
responds to the ideal condition of zero head loss, and that for m = 1 to the 
condition assumed by Favre of F = O. The unbalanced force (F) is directed 
downstream for m < 1 and upstream for m > 1. 

The experimental results, obtained by the writer and quoted by the author, 
indicate trends comparable to those for the curves of constant m. For 
D,/D = 1/4 and Dj/D = 1/2 the values of m are 1.38 and 0.95, respectively. 
The departure of the experimental curves from the theoretical curves for 


values of a less than 0.25 is due to the fact that for relatively small spines 


of Q) the type of mixing occurring results in a lower loss of head. For + 


= 1.0 the trend of the experimental curve is different from the theoretical 
curves, because the water entering the conduit through the lateral always ex- 
periences an increase in velocity; consequently, jet action does not occur and 
an entirely different flow pattern results. This trend is equally well supported 
by the Munich experiments. 
The energy equation gives a relationship for the loss of head between the 
lateral and the conduit as 
‘ 
4/29 (d) 
The experimental results correlating the change in _ for flow from the 
lateral to the conduit (n'y) with the discharge ratio 3 (shown in the author’s 


Fig. 8), can also be presented as in Fig. B. The results for the three cases, 
Le, D,/D equal to 1, 0.49, and 0.24, are presented. The two limiting 


Vv 2 
curves are defined by Eq. (d) with hg’ = 5 and hy = O. Most of the exper- 


2 


imental points lie between these two curves, particularly for values of + 


greater than 0.30 (which corresponds toz- = 0.3, 0.07, and 0.02 for 

D 

a = 1.0, 0.49, and 0.24, respectively). The points for the three cases have a 

tendency to form one continuous curve if plotted in this way. For higher 
An' 


values of aes , the values of ve og are very nearly zero. For values of 


Vv 
ma greater than 1.0, i is negative, a point discussed in some detail 
2 


by the author. 


The pressure measurements along the sides of the lateral just before the 
junction showed that the pressure on one side differed from that on the other 
side by very small amounts (approximately 0.005 ft. of mercury) in the case 
of the 1-inch lateral; this bears out the fact that for this case and F ~Q or 

= 1, These pressure differences were quite large (0.1 ft. of mercury) and 
opposite in sign for the 1/2-inch and 2-inch laterals, thus verifying the exist- 
ence ofa sizable unbalanced force F in the two directions for these two cases. 
From these observations and from the alternative plots (Figs. A and B), ad- 
ditional insight is gained in the interpretation and evaluation of the unbalanced 
force and head losses. 


RALPH W. POWELL,! M. ASCE—Despite the many papers that have been 
written on manifold flow, the subject has remained in a somewhat unsatisfac- 
tory state, and this treatment of the subject and the report on the results ob- 
tained at the Iowa Institute of Hydraulic Research is most welcome. The 
writer wishes however to call attention to one important mass of data on com- 
bining flow that had apparently not come to the author’s attention. 

In 1931 Naramoto and Kasai? reported on a series of 651 runs on combin- 
ing flow in cylindrical pipes with intersection angles of 15, 30, 45, 60, 70, 
75, 80, 85, and 90 degrees, and ratios of D, to D of 1.004, 0.815, 0.678, 
0.666, 0.581, and 0.503. Unfortunately for most of these runs only the total 
loss of energy is given, without telling how much of it was in the main pipe, 
and how much in the lateral. But for 29 of the runs for 90 intersections the 
complete data are printed, and from these the writer has prepared Tables 
l and 2, In general they show good agreement with the data reported in this 
paper, most of the points lying between the lines in Figs. 7 and 8 which rep- 
resent the Iowa and Munich experiments. 


TABLE 1 - Kyushu data with D,/D = 1.004 


2g Ah/Vv? 2 g hy /V? 2g hy / Vj 
0.135 0.450 - 0.480 - 26.7 
0.233 0.627 - 0.261 - 4.84 
0.341 0.939 0.100 0.82 
0.409 1,137 0.404 2.45 
0.484 1.161 0.368 1.59 
0.528 1.218 0.448 1.63 
0.564 1.268 0.510 1.64 
0.602 1.314 0.580 1.62 
0.653 1.352 0.648 1.54 
0.690 1.390 0.699 1.49 
0.704 1.396 0.684 1.43 
0.728 1.425 0.760 1.47 
0.745 1.447 0.808 1.48 
0.760 1.467 0.841 1.47 
0.767 1.453 0.846 1.4* 
0.774 1.461 0.878 1.46 
T. Prof. of Mechanics, The Ohio State Univ., Columbus, Ohio 
- “On the Loss of Energy at Impact of Two Confined Streams of Water” by 


Itaru Naramoto and Taijiro Kasai, Memoirs of the College of Engineering, 
Kyushu Imperial University, Vol. VI, No. 3, Fukuoka, Japan, 1931. 
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TABLE 2 - Kyushu data with D)/D = 0.503 


c 
3 


2g Ah/Vv? 2 g hy'/V? 2g hy'/ Vi 
0.135 0.396 - 0.184 - 0.642 
0.208 0.633 0.227 0.333 
0.270 0.797 0.656 0.574 
0.322 0.914 1.018 0.624 
0.368 1.043 1.442 0.680 
0.399 1.120 1.727 0.692 
0.435 1.193 2.085 0.703 
0.472 1.268 2.548 0.729 
0.509 1,132 2.990 0.734 
0.522 1.386 3.271 0.766 
0.567 1.496 3.813 0.756 
0.605 1.545 4.385 0.764 
0.598 1.543 4.294 0.765 


oon owt wow 


on 


JOHN K. VENNARD! A.M. ASCE, and DARRHL I. DENTONI,? J. M. ASCE— 
Congratulations to Professor McNown for his thorough analytical and experi- 
mental study of the manifold problem for simple basic geometries of the con- 
duit junctions. The writers have recently been involved in experimental work 
of a similar nature and are thus in a position to appreciate the excellent and 
concise presentation of theories and facts in a field which seems to have all 
too many “loose ends.” 

Since the completion of the thesis by James R. Barton [1], the senior 
writer has been intrigued by an apparently serious discrepancy between the 
results obtained at Iowa and at Munich for the case of dividing flow. This 
discrepancy is not evident on Professor McNown’s Fig. 4, in fact he makes 
the point that “the two series of experiments [i.e., lowa and Munich } indicate 
almost identical results.” From study of the Munich papers and comparison 
with Professor McNown’s interpretation of these on his Fig. 4, the writers be- 
lieve that the Munich curves have been mislabeled and the results are not only 
not “identical” but exhibit a trend opposite to those obtained at Iowa. This ap- 
parent discrepancy led to the junior writer undertaking an experimental pro- 
gram in an attempt to resolve the conflicting results obtained from these two 
carefully performed sets of experiments. Accordingly tests were conducted 
in the Stanford Hydraulic Laboratory on five square-edged outlets having axes 
at 90° to a main conduit of 1.276 inch diameter. Lateral sizes were selected 
which gave diameter ratios D,/D of 0.294, 0.392, 0.490, 0.642, and 0.830 
and the results plotted in the manner suggested by Professor McNown on his 
Figs. 4 and 5. Good agreement with the Iowa and Munich results was obtained 
in the case of Fig. 5 and fair agreement in that of Fig. 4. 


1, Professor of Fluid Mechanics, Stanford University, Stanford, California. 


2. Junior Engineer, Davis and Holman, Consulting Engineers, Stockton, 
California. 
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The Stanford tests were carried out in the sam, general manner as those 
at Iowa and Munich but differed in the one detail that no piezometer connec- 
tions were made in the laterals. The laterals consisted of smooth straight 
sections of pipe some 27 to 30 diameters in length discharging freely to the 
atmosphere. The diameters were measured (volumetrically) with extreme 
care and the hydraulic grade lines over the laterals established from the 
Blasius friction factor and the use of the generally accepted premise that the 
hydraulic grade line intersects the axis of the lateral in the plane of its outlet. 
The reliability of the Blasius friction factor, the correctness of the foregoing 
premise, and the adequacy of the above-mentioned length of lateral had been 
accurately confirmed by C. T. Chen in a series of tests connected with another 
experimental program. Many piezometer connections were used to establish 
the hydraulic grade line over the main conduit and the Blasius friction factor 
again used as a guide in selecting the best line through the tops of the piezom- 
eter columns. Here again the reliability of the Blasius friction factor was 
demonstrated by excellent agreement between a line of predetermined slope 
and points obtained experimentally. 

The Stanford results are presented in Figs. 1 and 2 of this discussion and, 
in view of Professor McNown’s discussion of his Figs. 4 and 5, are presumably 
self-explanatory. It is to be noted however that the trend of the curves con- 
firms that obtained at Iowa rather than that obtained at Munich and that the re- 
sults of Fig. 2 are in line with both the Iowa and Munich results. 

The lack of quantitative agreement between the results of the three test 
programs is not surprising since the quantities depend upon the relatively 
small vertical distances between extensions of hydraulic grade lines which are 
difficult to establish with great accuracy. It appears to the writers however 
that agreement is good enough and data now plentiful enough to develop a set 
of “standard values” for engineering use. Proceeding on the assumption that 
the accuracy of the results obtained at Munich, Iowa, and Stanford are about 
the same, “cross-plots” were prepared from points taken from the smooth 
curves obtained at the three institutions. This was accomplished by plotting 
4h/v?/2g against D)/D and hg'/V?/2g against D)/D for each even 0.1 of 
Q);/Q. These plots are not included here but the scatter of points on the first 
set was large and on the second relatively small. Because of the scatter it 
was judged that nothing better than straight lines were justified in the first set 
of crossplots; the best straight lines were fitted to the points by the method of 
least squares. For the second set of crossplots where scatter was less, curves 
were quite well defined and were drawn in by eye. The plots of Figs. 3 and 4 
were finally obtained; since these represent a composition and reconciliation 
of results from three separate test programs, they are offered as tentative 
“standard values” for engineering use. 

Although the writers do not wish to quibble over the symbols selected by 
Professor McNown, they feel that the use of the terms hy and h;' to symbolize 
local head losses is unfortunate since these symbols give the casual reader the 
impression that they embody some sort of “friction” phenomena. As an alter- 
native the writers would suggest that the upstream portion of the main conduit 
be designated by 1, the downstream portion by 2, and the lateral by 3. The 
term hg could then be replaced by h,, (or hy,,) and hg by h,, (for dividing 
flow) and h,, (for combining flow). This would seem to be more in line with 
generally accepted standard symbols. 
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CORRECTIONS TO THE DEVELOPMENT OF THE TURBULENT 
BOUNDARY LAYER ON STEEP SLOPES 
PROCEEDINGS-SEPARATE NO. 281 


The following figures are corrected versions of those which appear as 
Figures 7, 10, and 11 of Proceedings-Separate No. 281: 
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Fig. 7. - The Ratio £ as a Function of the Relative Roughness. 
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Fig. 10. - Graphical Solution of the Illustrative Problem. 
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Fig. 11. - The Shape Parameter as a Function of the 
Reynolds Number and Relative Roughness. 
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PROCEEDINGS-SEPARATES 
VOLUME 80 (1954) 


The technical papers published in the current calendar year are 
presented below. Technical division sponsorship is indicated by an 
abbreviation at the end of each Separate Number, the symbols refer- 
ring to: Air Transport (AT), City Planning (CP), Construction (CO), 
Engineering Mechanics (EM), Highway (HW), Hydraulics (HY), Irriga- 
tion and Drainage (IR), Power (PO), Sanitary Engineering (SA), Soil 
Mechanics and Foundations (SM), Structural (ST), Surveying and Map- 
ping (SU), and Waterways (WW) divisions. For titles and order cou- 
pons, refer to the appropriate issue of “Civil Engineering.” 


JANUARY: 379(SM),# 380(HY), 381(HY), 382(HY), 383(HY), 384(HyY),# 


385(SM), 386(SM), 387(EM), 388(SA), 389(SU),# 390(HY), 391(IR),# 
392(SA), 393(SU), 394(AT), 395(SA),2 396(EM), 397(ST). 


4Discussions grouped by Divisions. 
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